Double-dot charge qubit and transport via dissipative cotunneling 
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We investigate transport through an exotic "charge" qubit composed of two strongly capacitively 
coupled quantum dots (QDs), each being independently connected to a side gate which in general 
exhibits a fluctuating electrostatic field {i.e., Johnson- Nyquist noise). Two quantum phases are 
found: the "Kondo" phase where an orbital-Kondo entanglement emerges and a "local moment" 
phase in which the noise destroys the Kondo effect leaving the orbital spin unscreened and resulting 
in a clear suppression of the conductance. In the Kondo realm, the transfer of charge across the 
setting is accompanied by zero-point charge fluctuations in the two dissipative environments and 
then the I-V characteristics are governed by what we call "dissipative cotunneling" . 

PACS numbers: 73.23.-b, 72.15.Qm, 42.50.Lc 



The analogy between the Coulomb blockade effect in 
quantum dots (QDs) and the Kondo effect in heavy 
fermion systems has been emphasized for a decade 
Recently, fascinating quantum phase transitions were dis- 
covered in certain heavy fermion materials 0] , and a spin- 
boson-fermion impurity model has been proposed to cap- 
ture the quantum passage from local moment magnetism 
to strongly-correlated Kondo behavior Q . We seek to ex- 
plore the emergence of similar physics at the mesoscopic 
scale through a setup composed of a double-QD "charge" 
qubit immersed in a noisy electromagnetic environment. 
The central goal of this Letter is to focus on transport 
properties close to the quantum phase transition. We 
stress that a strong capacitive coupling between dots 
is essential for the fabrication of the double-QD charge 
qubit. Two capacitively coupled QDs have recently at- 
tracted much attention Q . One of the main motivations 
certainly lies in the fact that having coupled qubits is 
required for the controlled-NOT operation, one of the 
universal gates in quantum computing 5j. They are 
already known to exhibit non-trivial degeneracy points 
where quantum fluctuations can lead to exotic emergent 
states with orbital Kondo physics 0. Experimentally, 
coupled QDs can be built up, e.g., in semiconducting 
heterostructures and in nanotubes @- 

Model of two strongly capacitively coupled QDs subject 
to quantum noises. — The schematic circuit diagram of 
the setup under consideration is shown in Fig. 1. Each 
QD couples to a side-gate and each source of gate voltage 
is placed in series with an impedance Z\^.(uS) showing re- 
sistive behavior at low frequency u>. First, we will neglect 
the reservoirs of electrons (left and right leads) and seek 
to analyze the Coulomb blockade effect. At the Hartree- 
Fock level, the charging energy takes the form [3 

E ni n 2 = E cl [n g x - nif + E c2 (n g2 - n 2 f 

+E m (n g i - m) {n g2 - n 2 ) ; (1) 

we have E cl . 2 = e 2 C S2il /2(C sl C S2 - C 2 m ), E m = 
e 2 C m /(C^iC^ 2 -C m ), n gl . 2 = C gl , 2 V gl , 2 /e, and C S i,2 



FIG. 1: Left: Schematic setup of a double-QD charge qubit 
being connected in series to reservoir leads. The two QDs 
are capacitively coupled to side-gates and fluctuations of the 
gate voltages are included as a result of the finite resistances 
Ri.2 = Zi.2 (u> — 0) in the gate leads. Right: Honeycomb 
pattern of the charging energy for C m S> C g i where i = 1,2. 
Our concern is close to the degeneracy point n g i = n g 2 = 1/2. 



trons on each dot with respect to the equilibrium state 
at V g i,2 — 0. Our current interest lies in the degeneracy 



point n 



si 



lg2 



C g i_ 2 +C rn . In Eq. ffl. n\ and n 2 are the numbers of elec- 



1/2 and in the strong capacitive cou- 
pling limit C m > C g x,2, so that E i = E w < En = E 0() 
with En — Eqi — E m /2 ~ E cl ~ E c2 . The low-energy 
physics can thus be studied within the restricted Hilbert 
space in which only the (0,1) and (1,0) states are allowed 
(see Fig. 1). In this case, the double QD effectively re- 
duces to a charge qubit (two-level system) and one can 
introduce an orbital spin T z to describe the qubit, with 
the orbital spin up and down representing the (0,1) and 
(1,0) states respectively. It is convenient to introduce the 
operators Poi and P\q projecting on the states (0, 1) and 
(1,0) respectively. Through the fluctuation-dissipation 
theorem, the impedance Zi^ 2 (uj) introduces gate voltage 
fluctuations SV g i i2 . Since T z = (P i — Pio), the Hamil- 
tonian describing the charge fluctuations on the charge 
qubit due to 5V g i, 2 then can be written as ((SV g i(t)) = 0) 

H f i = e{5V gl -SV g2 )f z . (2) 

We plan to consider the case of Johnson-Nyquist noise 



2 



where the resulting gate voltage fluctuations are embod- 
ied by (SVgipSVgi^u = hRi,2 UJCOth(hui/2k B T); this 
can be generally attributed to the presence of ohmic re- 
sistances i?i t 2 in the gate leads. Below, we model the re- 
sistances i?i,2 in a microscopic fashion through two long 
dissipative transmission lines being composed of infinite 
collections of LiCu oscillators [Here % = 1,2 is the dot 
(and transmission line) index]. For a given gate lead, the 
resistive external circuit containing the transmission line 
has been previously introduced in Ref. I3. For the sake 
of clarity, here we follow Ref. by denoting the charge 
(fluctuation) operator on the capacitor C g i by Qo,i and 
the ones on the capacitors Cu between the inductances Li 
by Qi y i with I = 1, 2, and introducing the correspond- 
ing conjugate phase operators 4>i^. These transmission 
lines produce ohmic resistance Ri = y/ Li/Cu and are 
described by the bosonic Hamiltonian: 



Hp 



E 

i=l,2 



01 

2C a 



E 



14 h 2 (<f>i,i - 4>i-iaY 



2C'u 



+ 



2L, 



To diagonalize Hb analytically, we assume CV 



C a 



Using the transformation Qm = V2 dxcos[(l + 
l/2)irx]Qi(x) and = y/2 J Q dxcos[(l + l/2)nx]4>i(x), 
with [(f>j(x), Qk{y)] — ieS(x — y)6jk, and redefining Qi(x) 
and 4>i{x) in terms of boson operators a x %, we have 



(4) 



where u X i = 2ui C i sin(7rx/2) with uj C i — y/l/LiCu being 
a realistic high-frequency cutoff. The gate voltage fluc- 
tuations SVgi = Qoj/Cgi in Eq. (2) read ^1 



SVc 



a 1 



C, 



dx cos 



yJfiuj xl Cti (d xi + a+ ) . (5) 



It thus becomes clear that Eq. (2) introduces a coupling 
between the orbital spin and two boson baths. 

Now we incorporate electron tunneling between the 
two QDs. Throughout this paper we imply that a strong 
in-plane magnetic field has been applied, which allows us 
to restrict ourselves to spinless electrons. We still have to 
distinguish between large (metallic) QDs and small QDs. 
Large QDs are at the micron scale and then yield a very 
dense energy spectrum. In contrast, a small QD is at the 
nano scale and the energy level spacing is quite large. 

We first consider the case of large QDs. They are in 
the many-body state |$i$2) described by the free kinetic 
Hamiltonian H^ [ = J2i=i,2 J2 P , e Pidf Pi di Pi , where d lpi is 
the annihilation operator for an electron in the state pi 
on the ith dot. However, the Coulomb blockade effect 
analyzed above dictates that the tunneling can only take 
place between the (0,1) and (1,0) states. The key point 
is that one must use the projection operators Pqi an d 



P10 to actually project on | < I , i < I > 2) with charge number 
n\ = 0, ri2 = 1 and ri\ — l,ri2 — 0, respectively. The 
tunneling Hamiltonian then becomes 

^(large) = f £ rf + + ^ £ d+ p2 d lpi P 10 . (6) 



By explicitly introducing an auxiliary charge label \n1n2) 
to |$i ( I > 2) the Hamiltonian H t in Eq. © then turns into 



^■(large) _ ^ \ 



P1P2 



d t Pl d 2 P2 T- + E d t P2 d ipiT+- (7) 



P 1 P'2 



The tunneling matrix element t is assumed not to depend 
on the momenta p\ and p2- The T + operator - acting ex- 
clusively on \ri1n2) - ensures that each time an electron 
travels from dot 1 to dot 2 the charge state is explic- 
itly adjusted from 1 10) to |01> = f+|10). This causes a 
flip in the orbital spin which is accompanied by particle- 
hole excitations shared between the two dots. We find 
,(3) it appropriate to visualize as a transverse Kondo 

coupling along the lines of a single-electron box coupled 
to a reservoir lead The Kondo energy scale below 

which the transmission between dots becomes almost per- 



fect reads T^ rgc) = D e~ 1/(2|t|Aro) ; the ultraviolet energy 
cutoff is D ~ Pel, 2 and N is the density of states on 
the two dots which has been assumed to be symmetric. 

The total Hamiltonian for the two large dots in the 
absence of the leads can be summarized as 



H 



(lar| 
dot 



= H 



fl 



(large) 



H 



(o) 

dof 



(8) 



This is now clearly mapped onto a spin-boson- fermion 
model [l2^ . We find the following renormalization group 
(RG) equations: d\±/dl = X±\ z — gAj_/2, d\ z /dl — 
Aj_, and dg/dl = —g\\ with I = \h(Dq/D) denoting 
the scaling variable (D is the energy variable). Here, 
= 2|t|ATo, A 2 is the induced Ising-like coupling be- 
tween the orbital spin and the fermion degrees of free- 
dom, and g is an important physical parameter defined 
as g = 4 (i?i + R 2 ) /Rq (in the case of Cu = C g i)\ Rq = 
h/e 2 w 25.8fcf7 is the quantum of resistance. Like in Ref. 

0, we rewrite the RG equations in terms of the variables 
Aj. = 2|i|iVo and X z = \ z -g/2: initially, \ z = -g/2 < 0. 
Then, we recognize the well-known Kostcrlitz-Thouless 
(KT) equation flow d\±/dl = X±X Z and dX z /dl = A^, 
which reveals a quantum phase transition at X z = — Aj_, 

1. e., g = g c i = 4 1 i | TVo which experimentally corresponds 
to Ri + R2 = \t\NoRQ < Rq in the weak tunneling 
limit \t\No <C 1. For g < g c i, the charge qubit is still 
in the Kondo realm with the Kondo temperature being 
renormalized by noise: r£ argo) ( 5 ) « £) oe -i/(2|t|iv -g/2) . 
the boson baths have additive effects. For g > g c i, the 
noise effect becomes dominant and eliminates the orbital 
Kondo effect; the hopping between the two QDs becomes 
tiny and the current is strongly suppressed (see below). 

In contrast to the large QDs, the small QDs arc in a 
truly two-level state at the degeneracy point n g \ = n g 2 = 



3 



1/2 due to the large energy level spacing. Tunneling of 
electrons between the two QDs still induces an orbital 
spin flip however there is no particle-hole excitation in- 
volved. In consequence, the tunnel Hamiltonian between 
the two small dots takes the more conventional form \& 



H 



(small) 
dot 



H n +H B + (tT_ + h.c. 



(9) 



Such a model has been already investigated in great de- 
tail, including in the context of a small noisy dot em- 
bedded in a mesoscopic ring 01 ■ According to Ref. 0, 
this can be mapped exactly onto a transverse spin-boson 
model which in passing has been demonstrated to ap- 
ply to a wide range of two-level state problems including 
a superconducting qubit 0] and a quantum dot con- 
nected to Bose- Einstein condensate reservoirs 15]. The 
phase diagram can also be derived from an analogy with 
the single-impurity (anisotropic) Kondo model |13| . This 
already ensures some universality in the transport fea- 
tures independently of the precise sizes of the QDs. Here 
the Kondo energy scale T ^ mall \g) = |i|{|i|/A)} 9/(1 ~ s) 
only decays as a power-law with g and that the KT phase 
transition should arise in the vicinity of g = g cs ~ 1 |l0| . 

Now, we are in a favorable position to accurately study 
transport properties through the double-QD charge qubit 
in the quantum limit of T = 0. We treat the case of the 
large dots and the small dots in order. 

Transport through large QDs. — We connect the two 
QDs in series to two leads (source and drain) and add a 
voltage drop V across the two leads as shown in Fig. 1. In 
the restricted Hilbert space, at small energy eV, single- 
electron sequential tunneling through the dots is expo- 
nentially suppressed 0] . The leading order contribution 
to the tunneling here comes from the cotunneling process: 

4^°) = t > J2 E (4 k2 c lkl d+ P1 d2 P2 f- + h.c.) , (io) 

kik 2 P1P2 

where annihilates a spinless electron in the state fcj on 
the left (i = 1), right (i = 2) lead; t' ~ t L t R /[Eu-Eoi} = 
2tL,tn/E m with ti, (tn) being the tunneling matrix el- 
ement between the left lead and dot 1 (the right lead 
and dot 2). Since H^ 1 ^ involves four fermion opera- 
tors, this is irrelevant in the RG spirit and can always be 
treated perturbatively. 

To compute the tunneling current, it is convenient to 
absorb the coupling to the quantum noise into the tun- 
neling term through the canonical transformation hi = 
exp{iT z (<j) 0A - 4> ,2)}, where (j> ,i = f /* dt' SV gi (t')dt' . 
The Hamiltonian (|10fl then turns into: 



H, 



(large) 



= t'A + h. 



^=EE c ik 2 c ikidt Pl d2 P2 T- 

kik 2 P i P2 

By using linear response theory, the tunneling current 
then can be expressed as 0] 



(11) 

0,1-00,2) (\2) 



where (• • ^ret means the retarded correlation function 
computed with respect to the states of H^ s °\ It is 
important to note that the Hamiltonian in (|1C)|> alone 
does not drive a current because it conserves the total 
number of electrons in the left lead and in the left dot 
0. It is therefore necessary to have the tunneling term 



between the two QDs, H. 



(large) 



in order to obtain a non- 



vanishing tunneling current. We immediately infer that 
in the boson phase g > g c i, the current is negligible as 
a consequence of the suppressed tunneling between the 
two dots due to the fluctuating electromagnetic fields. 
Now we focus on the Kondo regime g < <? c ;. As dis- 



cussed before, the Kondo coupling X±(l) in H 



(large) 



flows 



to the strong coupling limit resulting in a perfect trans- 
mission between the two dots. To the leading order, we 
then expect (E PlP2 d f pi d 2p2 f„ ■ ~£ P3Pi d + Pi d lpa f + ){uj) ~ 
5(ui) as a result of the formation of an Abrikosov-Suhl 
resonance at zero frequency (lu = 0) which is the signa- 
ture of a Kondo ground state, and thus 



I L = 2ne\t'\'N L N R / ^ / dui / dcj 2 /(£) 

J~oo JO JO 

x [1 - /(£ + eV + hut + ftw 2 )FiK) Vzfa). (14) 

Here /(£) is the Fermi function (At low eV, we can re- 
strict ourselves to equilibrium states HI Ea), N L and 
Nr are the density of states in the left and right leads, 
and we introduce the following Fourier transform 



dte iwt( e i4> Oti (t) e -4 o>i (0)} 



[exp(-a l7 o)/r(a t )] ^"V^S 



(15) 



I L = -2e\t'\ 2 Im(A-A + ) Iet (-eV), 



(13) 



for lj < u c , where a, = 2RJRq (so cti + a 2 = g/2), 70 ~ 
0.5772 is Euler's constant, and T(x) the Gamma function. 
For convenience, we have assumed that lu Ci — uj C2 = lu c . 
The quantity Tjiujj) has a well-defined mesoscopic inter- 
pretation [lfllla,liJi : this represents the probability that 
an electron which tunnels onto the dot i also creates a 
charge excitation with an energy fkui in the dissipative 
environment i. From Eq. (|14fl it becomes apparent that 
the cotunneling of electrons across the double dot must 
be mediated through excitations in the two environmen- 
tal bosonic modes; this will be referred to as "dissipative 
cotunneling" . Inserting Eq. (|15|l into Eq. (|14Jl yields 

I L = c L [exp(-<7 7 o/2)/r(2 + g/2)] (eV/hu c ) 1+9/2 (16) 

at T = 0, where cl = \t'\ 2 NLNpfujj c /eRQ. Again, we 
see that the two noises have additive effects. In Fig. 2, 
we plot II normalized to cl as a function of eV/fwj c at 
various g. We discern that in the Kondo regime the I— V 
characteristic curves of the double-QD charge qubit can 
be viewed as an interesting generalization of the ones of 
the noisy tunnel junction 1 1 St Il9| and QD ^^|- The de- 
creasing current with growing g at small voltage can be 



4 




eV/(0 

c 

FIG. 2: Normalized currents Il/cl and Is/cs as a func- 
tion of the dimensionless bias voltage eV/u c (h = 1) in the 
Kondo regime. Curves from top to bottom correspond to 
g/2 — ot\ + Q2 = 0, 0.1, 0.2, 0.5. The inset is a log-log plot of 
the conductance for large dots (solid line) and for small dots 
(dashed line) versus \t'\ in the Kondo regime in the specific sit- 
uation where g = and Tkl > T^ mall) (g) Qt\/D = 0.0001). 



associated with the more and more suppressed noise spec- 
trum at small u. In the absence of the noises, the con- 
ductance obeys Gl — W\ 2 Nt,NrI Rn similar to a large 
QD at the Coulomb blockade peaks |l(| . 

Transport through small QDs. — When the two small 
QDs are connected in series to the leads, the cotunneling 



Hamiltonian becomes H, 

.-(small) 



(small) 



*' £ felfe2 (cL 2 c ifci T - 



h.c). Here Hf ' can also generate a Kondo effect at 
the energy scale Tkl = -Doe -1 /' 2 '* ' Nl ^ (we have assumed 
Nl = Nr — Ni). For quite reasonable t, in principle, we 
should have r p^ mall "> (g) ^> r p R l w hich guarantees the va- 
lidity of our perturbative treatment of _ff^ small \ i.e., Eq. 
(fTCfl found previously for the large dots is then applica- 
ble. Eventually, by converging to the extreme limit of 
very small t, H^ mai1 ' could dominate the orbital Kondo 
physics and then the leading order term in the current 
becomes quadratic in t || . The current at T = reads 



Is = c s [exp(- 97 o/2)/r(2 + g/2)](eV/Tuj c 



,1+9/2 



(17) 



where cs — \t/TKL\ 2 nLo c / eRQ. For a = 0, we recover 
the conductance G s = \t/T KL \ 2 /R Q Q. Both I L in Eq. 
(|16|l and Is in Eq. (|17fl have the same power law, albeit 
due to different orbital Kondo physics. A comparison be- 
tween Gl and Gs versus \t'\ has been performed in Fig. 2 
(inset). Note that orbital Kondo physics differs from the 
conventional spin Kondo effect in which the perfect trans- 
mission usually leads to the conductance quantum 2/ Rq 
(the factor 2 is attributed to the two spin channels). 

In conclusion, we have carefully investigated the ef- 
fect of gate voltage fluctuations on transport through a 



double-QD charge qubit. We have obtained striking I-V 
characteristics as a result of the prolific combination be- 
tween orbital Kondo entanglement and dissipative cotun- 
neling (the transfer of electrons through the structure has 
to be conducted by bosonic excitations in the two electro- 
magnetic environments). This generalizes previous noise 
studies on the single tunnel junction and QD. When the 
coupling to the dissipative environments becomes large 
enough the orbital Kondo effect is suppressed and the 
orbital spin acts like a free moment which completely de- 
pletes the current (at low V). A similar analysis can be 
performed with superconducting leads and dots. 
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